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' Abstract 

| We discuss the advantages of using the approximate quantum Fourier trans- 

■ form (AQFT) in algorithms which involve periodicity estimations. We analyse 

""^5 \ quantum networks performing AQFT in the presence of decoherence and show 

that extensive approximations can be made before the accuracy of AQFT (as 
compared with regular quantum Fourier transform) is compromised. We show 
that for some computations an approximation may imply a better performance. 
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u • 1 Introduction 



In the course of history many ingenious mechanical, acoustic and optical devices have 
been invented for performing Fourier transforms |ffl] (including nature's own such as 
the human ear) . Most of them are now of merely historical interest since the arrival of 
the computer-based algorithm known as the fast Fourier transform (FFT) |2|, |J| which 
efficiently computes the discrete Fourier transform. The FFT algorithm can also be 
phrased in terms of quantum dynamics, i.e., in terms of unitary operations performed 
by a quantum computer on quantum registers. Indeed, all known quantum algorithms 
employ the quantum version of Fourier transforms, either explicitly or indirectly. It 
is used for the periodicity estimation in the Shor algorithm Q and its approximate 
version (the Hadamard transform) is commonly used to prepare quantum registers 
in coherent superpositions of different values. 

In this paper we analyse the performance of the quantum Fourier transform 
(QFT) in the presence of decoherence. In particular we show that as far as the peri- 
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odicity estimation is concerned the approximate quantum Fourier transform (AQFT) 
can yield better results than the full Fourier transform. 

In the following we use some terms which were originally adopted from the 
classical theory of information and computer science and became standard in the lore 
of quantum computation. More detailed descriptions can be found e.g. in Refs. [||, o, 

01 and in some recent reviews ||. 

• A qubit is a two-state quantum system; it has a chosen 'computational basis' 
{|0), |1)} corresponding to the classical bit values and 1. Boolean operations 
which map sequences of O's and l's into another sequences of O's and l's are 
defined with respect to this computational basis. A collection of L qubits is 
called a register of size L. 

• Information is stored in the registers in binary form. For example, number 6 is 
represented by a register in state |1) ® |1) <8> |0). In more compact notation: \a) 
stands for the direct product <8> \cll~2) . . . |oi) ® |oo) which represents a 
quantum register prepared with the value a = 2°ao + 2 l a\ + . . . 2 L ^ l a^\. 

• A quantum logic gate is an elementary quantum computing device which per- 
forms a fixed unitary operation on selected qubits in a fixed period of time. 

• A quantum network is a quantum computing device consisting of quantum 
logic gates whose computational steps are synchronised in time. The outputs 
of some of the gates are connected by wires to the inputs of others. The size 
of the network is its number of gates. 

• A quantum computer will be viewed here as a quantum network (or a family of 
quantum networks). Quantum computation is defined as a unitary evolution of 
the network which takes its initial state "input" into some final state "output" . 

Our presentation starts with a brief mathematical introduction to the approxi- 
mate discrete Fourier transform which is followed by the description of its quantum 
implementation in terms of quantum networks. Then we analyse how the perfor- 
mance of the QFT in the periodicity estimation is affected by the approximations in 
the algorithms and by decoherence. We also comment on possible simplifications in 
practical implementations of quantum networks effecting the QFT and AQFT. The 
quantum algorithm for the fast Fourier transform which we use in this paper was 
originally proposed by Coppersmith and by Deutsch (independently) ||. 

2 Discrete Fourier Transforms 

The discrete Fourier transform is a unitary transformation of a s-dimensional vector 



{/(0), /(l), /(2), . . . , f(s - 1)} defined by: 
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where /(a) and /(c) are in general complex numbers. It can also be represented as 
a unitary matrix 
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where u> = exp(27ri/s) is the sth root of unity. In the following we assume that s is 
a power of 2, i.e., s = 2 L for some L; this is a natural choice when binary coding 
is used. The approximate discrete Fourier transform can be conveniently described 
when we write the product ac in the exponent on the r.h.s. of Eq. ([!]) in the binary 
notation. Writing 
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we obtain 



ac 



a c + (a ci + aic )2 + (a c 2 + aici + a 2 c )2 2 + . . . 



+ (a c L _i + ... + a L _ic )2 L - 1 +C>(2 



(4) 



Because lo x = 1 for x > s, the terms 0{2 L ) do not contribute to the transform, and 
the term exp(27riac/2 L ) in Eq. (0) can be expressed as 

exp(27rmc/2 L ) = exp(27rz(a c )/2 L ) exp(27rz(a ci + aic )/2 L_1 ) . . . 

x exp(27rz(a CL_i + . . . + a L _ic )/2). (5) 

Beginning from the right of this expression, the arguments in the exponentials become 
smaller and smaller. In the approximate Fourier transform parameterised by an 
integer m, the L — m smallest terms are neglected. In all the remaining terms the 
arguments are then multiples of 27r/2 m . The 2 m th root of unity becomes the basic 
element of the approximate Fourier transform as opposed to 2 L th root of unity which 
used in the ordinary Fourier transform. (The ordinary Fourier transform is obtained 
for m = L; when m = 1 we obtain the Hadamard transform, for which all terms but 
the last one are dropped.) 

The quantum version of the discrete Fourier transform is a unitary transfor- 
mation which can be written in a chosen computational basis {|0), |1), . . . , \ s — 1)} 
as 



QFT S 
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s-l 



—= ^ exp(2-7riac/.s) |c) 
v s c =o 



(6) 



More generally, QFT S effects the discrete Fourier transform of the input amplitudes. 
If 

QFT s :£/(a)|a)^£/(c)|c), (7) 

a c 

then the coefficients /(c) are the discrete Fourier transforms of /(a)'s. This defini- 
tion can be trivially extended to cover the approximate quantum Fourier transform 
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(AQFT). We will analyse the approximations involved in the AQFT in terms of 
computational networks. 



3 Quantum Networks for AQFT 



Quantum networks for AQFT can be constructed following the description of the 
fast Fourier transform algorithm (as described by Knuth [jlO| ) . This efficient classical 
algorithm needs to be re-expressed in terms of unitary operations ||. The con- 
struction requires only two basic unitary operations. The first operation is a one-bit 
transformation A4 (one-bit gate) that acts on a qubit qi of the register and effects 
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The diagram on the right provides a schematic representation of the gate acting on 
a qubit q. The second operation is a two-bit gate Bjk that effects 
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where 9jk depends on the qubits j and k on which the gate acts and equals Ojk = 
Tr/2 k ~i . The transformation Bjk is an elementary two qubit operation which affects 
only states with a 1 in both position j and k regardless the state of the remaining 
qubits. 

The QFT on a register of size 1 reduces to a single operation A performed on 
a single qubit (cf. equation (^) for s = 2). The extension of the QFT network to a 
register of any size L follows a simple pattern of gates A and B which can be seen 
in Fig. [lj. It shows the QFT network operating on four qubits which can be written 
as the sequence of the following 10 elementary operations (read from left to right) 



{A3 ) (-B23 A 2 ) {Bi 3 B12 A\ ) (B 3 B 2 Bqi A ) . 



(10) 



The bit values at the output should be read in the reversed order (see Fig. p]). 

The number of gates needed to complete the QFT grows only as a quadratic 
function of the size of the register: for a transformation on a L qubit register, we 
require L operations A and L{L — l)/2 operations B, in total L{L + l)/2 elementary 
operations. Thus the quantum QFT can be performed efficiently. 

The AQFT of degree m is represented by a similar network in which the two- 
bit gates that act on qubits which are far apart (in the register) are neglected, i.e., 
those operations Bjk for which the phase shift 6jk = 7r/2 fc_ -' < 7r/2 m for some m 
such that 1 < m < L are dropped (cf. Eq.@). In that case, we need L operations 
A, and (2L — m)(m — l)/2 operations B, which is an improvement on the QFT case 
since m < L. In Fig. ^ we show the m = 2 AQFT network counterpart to the QFT 
network shown in Fig. ||. The matrix elements of the QFT and the AQFT differ by 
a multiplicative factors of the form exp(ie) with |e| < 2irL/2 m . The execution time 
of the AQFT grows as ~ Lm. 
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4 Estimating Periodicity 



The quantum Fourier transform, like the ordinary Fourier transform, is a powerful 
tool for uncovering periodicities. Any periodicity in probability amplitudes describing 
a quantum state of a register in a computational basis can be estimated (with some 
probability of success) by performing the QFT followed by a measurement of the 
register in the computational basis. The result is obtained by reading the qubits of 
the register in the reversed order. 

For example, an interesting periodic state which plays an important role in 
Shor's quantum factoring algorithm can be written as 

2^ — 1 

V-A/ a=0 

where M is an appropriate normalisation factor and 

f(a) = Si, a mod r- (12) 

It is the state of a quantum register of size L in which the probability amplitudes 
/(a) occur with periodicity r and offset I. If this offset is unknown, a measurement 
performed in the computational basis cannot reveal r or any of its integer multiples 
directly. This is illustrated in Fig. §(a). However, if we perform the QFT on the 
register first and subsequently measure its state we obtain number c which, with 
probability greater than 4/-7T 2 , is a multiple of 2 jr regardless the offset I (cf. Ap- 
pendix). The probability is not equal to unity, because the finite size of the register 
leads to a "broadening" of the Fourier-transformed data, as illustrated in Fig. [|(b) . 
(This is because 2 L jr is not necessarily an integer, and the quantum register can 
have only integer values; this is discussed in detail in the Appendix). 

For the AQFT the corresponding probability, in the limit of large L, satisfies 

Prob^sm 2 ^-]. (13) 

This result is derived in the Appendix. The effect of the approximation is illustrated 
in Fig. |D where we plot the modulus of the amplitude of the transformed state 
(with I = 9 and r = 10) for the AQFT of different orders m. Fig. || shows how the 
phase of the transformed state becomes corrupted when m becomes smaller. 

If the quantum Fourier transform forms a part of a randomised algorithm then 
the computation can be repeated several times in order to amplify the probability of 
the correct result. In such cases the performance of the AQFT is only polynomially 
less efficient than that of the QFT. For example, consider Shor's quantum factoring 
algorithm and substitute the AQFT for the QFT. In order to obtain a correct factor 
with a prescribed probability of success, we have to repeat the computation several 
times. Let k and k' be the number of runs respectively with the QFT and the AQFT 
so that we obtain the same probability of getting at least one correct result, i.e., 

l-(l-p) k = l-(l- p 'f. (14) 
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Here p = 4/ir 2 and p' = ^ sin 2 are the corresponding probabilities of success 

in a single run. The ratio k' jk scales as 




for some C (the upper bound is found graphically). This shows that in the quantum 
factoring algorithm the AQFT is not less efficient than the ordinary QFT, i.e., the 
ratio k'/k scales only polynomially with L/m. Moreover, we will show that in the 
presence of decoherence the AQFT can perform better than the QFT even in a single 
computational run! 



5 Decoherence 

Quantum computation requires a controlled, quantum-mechanically coherent evo- 
lution at the level of individual quantum systems such as atoms or photons. This 
imposes severe requirements on quantum computer hardware. If we are to harness 
the unique power of quantum computers, such systems will have to be manufactured 
with unprecedented tolerances and shielded from noise to an unprecedented degree. 
Even a minute interaction with the environment will lead to a non-unitary evolution 
of the computer and its state will, in general, evolve into a mixed state described 
by a reduced density operator p, which is obtained from the density operator ptotai 
of the total computer+environment system by taking a trace over all the quantum 
states of the environment: 

P — Tr cnv i ronmen t (ptotal)' (16) 

Consider, for example, a quantum register of size L which is prepared initially in some 
pure state and then left on its own. As time goes by, the qubits become entangled 
with the environment. Both the diagonal and the off-diagonal elements of the density 
matrix p (expressed in a computational basis) are usually affected by this process 
(cf. |ll], 0|). The rate of change of the diagonal and the off-diagonal elements 
depends on the type of coupling to the environment, however, there are realistic 
cases where the dissapearance of the off-diagonal elements, known as decoherence, 
takes place on much faster time scale. In this case a simple mathematical model of 



decoherence has been proposed [13|. It assumes that the environment effectively acts 
as a measuring apparatus, i.e., a single qubit in state co|0) + ci|l) evolves together 
with the environment as 

(co|0)+ci|l))|o) -^co|0)|a ) + ci|l)|ai), (17) 

where states \a), |ao), |ai) are the states of the environment and |ao), \a±) are usually 
not orthogonal ((ao|ai) ^ 0). The elements of the density matrix evolve as 

Pij {0) = q(0)c*(0) — ► pij{t) = c J (t)c*(t)(a J (t)|a J (t)), i,j = 0, 1. (18) 

States |do) an d \a\) become more and more orthogonal to each other whilst the 
coefficients {cj} remain unchanged. Consequently the off-diagonal elements of p 
disappear due to the (ao(t)\ax(t)) factor and the diagonal elements are not affected. 
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There is an alternative way of thinking about this process. The environment 
is regarded as a bosonic heat bath, which introduces phase fluctuations to the qubit 
states, i.e., it induces random phase fluctuations in the coefficients cq and c% such 
that 



c | 0) + ci | 1) 



c e 



-i(f> | 



0) + cie^ | 1 



The direction and the magnitude of each phase fluctuation 
following the Gaussian distribution 



P 



1 



2tt8 



exp 
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(19) 

is chosen randomly 
(20) 



where the distribution width 5 defines the strength of the coupling to the quantum 
states of the environment. The elements of the density matrix p are then recon- 
structed as pij = (qc*), where the average () is taken over different realisations of 
the phase fluctuations within a given period of time (cf . |l2| ) . The diagonal elements 
do not depend on (j), whereas the off-diagonal term (cocie* 2 ^) averages to zero for a 
sufficiently long period of time. 

The latter approach to decoherence is very convenient for numerical simulations 
and was chosen for the purpose of this paper. It is similar to the Monte Carlo wave 



function method used in quantum optics [14 



6 AQFT and Decoherence 

We have analysed decoherence in the AQFT networks assuming that the environment 
introduces a random phase fluctuation in a qubit probability amplitudes each time the 
qubit is affected by gate B. In our model we have not attached any decoherence effects 
to gate A. In most of the suggested physical realisations the single qubit operations 
are quite fast, whereas the conditional logic needed in two-qubit operations is often 
much harder to produce, which makes these operations slower than the single qubit 
ones, and often much more susceptible to decoherence. For instance, in the ion trap 
model proposed by Cirac and Zoller fl5| single qubit operations require only one laser 
pulse interacting with one atom, whereas in a two-qubit operation two subsequent 
laser pulses are needed, and the atoms involved must form an entangled state with a 
trap phonon mode between the pulses. In any case, introducing decoherence to gate 
A operations and 'wires' in the network would not affect our results much, only the 
time scale for decoherence would change. 

We have quantified the performance of the AQFT by introducing the quality 
factor Q. It is simply the probability of obtaining an integer which is closest to 
any integer multiple of 2 L /r, when the state of the register is measured after the 
transformation. In the decoherence free environment analysed in Sec. || we obtain 
Q = 1 for integer values of 2 L /r and for a randomly selected r the quality factor Q 
is of the order 4/vr 2 for the QFT and of the order of ^ sin 2 \ for the AQFT of 
degree m. 

In Fig. ^ we show how the quality factor Q behaves as a function of m and 5 
(which characterises the strength of the coupling to the environment). For S > 
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the maximum of Q is obtained for m < L. Thus in the presence of decoherence one 
should use the AQFT rather than the QFT. 

This 'less is more' result can be easily understood. The AQFT means less 
gates in the network and because each B gate introduces phase fluctuations the 
approximate network generates less decoherence as compared to the regular QFT 
network. By decreasing m we effectively decrease the impact of decoherence. On the 
other hand decreasing m implies approximations which reduces the quality factor. 
This trade-off between the two phenomena results in the maximum value Q for 
m G [l,L]. 

It is worth pointing out that for 5 = (no decoherence) Q remains almost 
constant for those values of m that satisfy the lower bound condition (derived in the 
Appendix) 

m>log 2 L + 2. (21) 

and when 5 > the optimum m is found near this lower bound. In Fig. ^ we also 
show how Q decreases rapidly with L in the QFT network (although there is not 
enough data in the figure to determine if it really decreases exponentially). 

Our simulations were performed for ensembles which consisted typically of a 
one to two thousands individual realisations. 



7 Conclusions and comments 

We have analysed the approximate quantum Fourier transform in the presence of de- 
coherence and found that the approximation does not to imply a worse performance. 
On the contrary, using the periodicity estimation as an example of a computational 
task, we have shown that for some algorithms the approximation may actually imply 
a better performance. 

Needless to say, there is room for further simplifications of the quantum Fourier 
transform which may lead to at least partial suppressing of unwelcome effects of 
decoherence. For example, if the QFT is followed by a bit by bit measurement of the 
register then the conditional dynamics in the network can be converted to a sequence 
of conditional bit by bit measurements (cf. |16[j). Here, we wanted to show that there 
are cases where quantum networks that are composed of imprecise components can 
guarantee a 'pretty good performance'. This topic has also a more general context; 
it has been shown that reliable classical networks can be assembled from unreliable 



components [17|. It is an open question whether a similar result holds for quantum 
networks. 
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Consider the quantum state 



I*) = —rrf E /(«)!«)' ( 22 ) 

where 

/(«) = Si, a mod r- (23) 

Here /(a) is a periodic function with period r <C 2 L and offset Z, which is an arbitrary 
positive integer smaller than r; see Fig. ||(a). The normalization factor is equal to 
the number of non-zero values /(a): M = [2 L /r]. Because r -C 2 L we use from 
now on [2 L /r] ~ 2 L /r. The state (22) plays an important role in the Shor quantum 



factorisation algorithm (the algorithm enables us to factorise an integer N by finding 
r such that x r = (1 mod N) for some x coprime with N — r is estimated from a 
quantum computation that generates a state of the form (P2|).) 
Applying QFT to this state we obtain 



l^>=E/( c )l c >> (24) 



where 



It 



Kc) = \r E exp(2vri(jr + /)c/2 L ). (25) 
1 3=0 

The probability of seeing an integer c is then 

2 

(26) 



Prob(c) = \f(c)\ 2 = ^ 



2 L /r-l 

E exp(2vrij(rc mod 2 L )/2 i 
3=0 



As it can be seen from Fig. |^(b) the peaks of the power spectrum of f(a) are centered 
at integers c which are the closest approximation to multiplies of 2 L /r. 

Let us now evaluate Prob(c) for c being the closest integer to X2 L /r, i.e., c = 
[X2 L /r\. By definition c must satisfy 

1 2 L 1 

-i< 5 -A^<i. (27) 

We define 6c = 27r(rcmod 2 L ) so that Prob(c) now involves a geometric series 
with ratio exp(i# g ). By viewing these terms as vectors on an Argand diagram it 
is clear that the total distance from the origin decreases as 9c increases. Hence 
Prob(c) > Prob( c with largest allowed 5 ). Let us denote by # max the largest allowed 
9c- By Eq. (p7|), 9 max < irr/2 L and summing the geometric series with # max = vrr/2 L 
(see Fig. ||) we obtain 

r 1 4 1 

Prob(c) > t ^ ~ (28) 

2 sin 2 vr 2 r 



(i 



2^ 
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where we have used the fact that r/2 L is small. Since there are r such values c, the 
total probability of seeing one of them is 



Prob > 4/tt 2 . 



(29) 



By performing this measurement several times on different states \^f) (each one with 
possibly different /), one gets with high probability values Co, ci,... that are the 
closest integers to uq2 l /r, n\l L /r, . . . and which allow to calculate r [cf. inset in 
Fig. 1(b)]. 

We estimate now the probability of measuring one of the desired values c when 
the AQFT of order m has been performed instead of the QFT. The difference between 
the QFT and the AQFT of order m is in the arguments of the exponentials in Eq. 
(p6[). The phase difference for each term in the sum is 



L-l 



A(a,c) = -r 



ac 



a j c k 2 j ~ 



(30) 



where 

£ = {(j, k) | < j, k < L - 1, L - m < j + k < L - 1} . (31) 

The probability to measure |c), where c is the closest integer to one of the r 
values n2 L /r now becomes 



Prob A (c) = 



2 L /r-l 

Y exp(2nij(rc mod 2 L )/2 L - iA(jr,c)) 



(32) 



This is the same summation as is involved in the QFT, except that in the case of the 
AQFT, each vector of the Argand diagram of Fig. B(a) may be rotated by an angle 



A(jr, c), as shown in Fig. 
Vi, A (a, c) is equal to 



In the worst case, when a 



1, i.e., en 



2vr 



m 



l + 2 r 



1 



(33) 



However, for any other values of a and c, < A(a, c) < A max . 

We are interested in the lower bound for the probability so we assume that the 
vectors in the Argand diagram fill one half of the circle (# max = vrr/2 L ) as illustrated 
in Fig. |§|(b). The approximation allows to rotate each vector by the maximum angle 
A max - The minimum of the probability is obtained when half of the vectors are 
rotated by A max , see Fig. ^(b). In this case vectors in two areas of size A max cancel 
each other, and all we have to do is to calculate geometrical sums of the vectors in 
the two areas of size tt/2 — A max . In an area of that size there are (j} — A ™ ax ) 
vectors, since the total number of vectors is 2 L jr. Note that because 2 L r, we can 
assume that [2 L /r] ± 1 ~ 2 L jr. The square of the geometric sum then becomes 



' 1 Amax * 

,2 IT 

E 



-i 



ex P 1 ijrJ 



sin 



(2(2 Amax) 

(fir) 



(34) 
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The two sum vectors in the two areas of size tt/2 — A max are of equal length 
and orthogonal to each other, so the square of their sum vector, contributing to the 
total probability, is twice the value given by (pi)- Finally we obtain 



r 2 sin 2 (± (| - A max ) J g „ /l /vr \\ 
Pr°b A >2^ ^(^) ^ -,sij(-(--A mm )). (35) 

For 

A max — this expression reduces to the result derived for the QFT: 

2 4 



Prob A >- I sin^ I j=- f , (36) 

and for A max = tt/2 we have Prob^ > 0. To avoid Bb zero probability, A max must 
always be bounded by 

A max = ^(L-m-l + 2 m - L ) < |, (37) 
which for large L implies 

m > log 2 L + 2. (38) 

Eq. (|2l| ) gives a lower bound to the order of the AQFT performed on a register 
of length L, if one wants to have a non-zero probability of success in measuring a 
value c. Simple geometric considerations also show that A max < tt/2 is a limit for a 
non-negligible probability: for A max > tt/2 the vectors in the Argand diagram can 



be rotated so that there is no longer any constructive interference, see Fig. 1C. 
For large L we can write 

2ir 

A max ~— (L-m). (39) 
If we use the lower bound for m (^), we obtain 

A max < - (1 - -J , (40) 
which allows to write the probability ( |35| ) in a simple form 

Prob A >^sin 2 ( IT j. (41) 
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Figure 1: QFT network operating on a four-bit register. The phases Ojk that appear 
in the operations Bjk are related to the "distance" of the qubits (j — k) and are given 
by 9jk = 7r/2- y_fc . The network should be read form the left to the right: first the 
gate A is effected on the qubit 0,3, then B^^) on 0-2 and 0,3, and so on. 
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Figure 2: AQFT network operating on a four-bit register when m = 2. 



Figure 3: (a) Function /(a) in Eq. (12). The parameters are / = 8, r = 10 and the 
number of bits in the calculation is L = 9 so that numbers up to 2 L -1 = 511 can be 
encoded, (b) |/(c)|, obtained from f(a) by a QFT. The inset shows C3 = 155 which 
is the closest integer to 3 • 512/10 = 153.6. (See the Appendix for more details). 

Figure 4: Different orders of approximation in the AQFT performed on a state 
for which f(a) = 6 9) a mod 10 . 

Figure 5: As Fig. ||, but showing the phase of the amplitudes, i.e., arg(/(c)). 

Figure 6: The quality factor Q as a function of m for selected values of 8. The register 
sizes are (a) L = 9 and (b) L = 16. Statistical errors are too small to be represented 
on the graph. 
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Figure 7: The quality factor Q as a function of the register size L for QFT, with 
varying levels of decoherence, from 5 = 0.1 (top line) to 5 = 0.5 (bottom line). 
Statistical errors are too small to be represented on the graph. 



Figure 8: (a) Argand diagram corresponding to the sum of the phases that appear 
in the expression of f(c) for c close to one of the values n2 L /r. Prob(c) is the norm 
of the vector resulting from the sum of each vector in the diagram, (b) Prob(c) is 
bounded by the worst case situation in which we have taken # max instead of 6 5 , in 
this case the phases lie on an interval [0, it] on the Argand diagram and a closed form 
expression can be found. 



Figure 9: (a) Argand diagram for the AQFT. Vectors are rotated by an angle 
A(jr, c). (b) To obtain a closed form for a bound for Probyi(c) we consider the worst 
case in which half of the phases pick up a factor A max . 



Figure 10: In the case that the order m of the AQFT is such that A max > it/2, 
the individual phases can get scrambled in such a way that there is no constructive 
interference effect. The probability Probyi(c) can become vanishingly small and the 
AQFT of order m is inefficient. 
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